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Abstract 

By considering the master equation of the partially asymmetric diffusion 
process on a one- dimensional lattice, the most general boundary condi- 
tion (i.e. interactions) for the multi-species reaction-diffusion processes 
is considered. Resulting system has various interactions including diffu- 
sion to left and right, two-particle interactions A a Ap — > A~,A& and the 
extended n-particle drop-push interactions to left and right. We obtain 
three distinct new models. The conditions on reaction rates to ensure the 
solvability of the resulting models are obtained. The two-particle condi- 
tional probabilities are calculated exactly. 



1 Introduction 

The understanding of non-equilibrium statistical physics is still much more in- 
complete than that of equilibrium theory, due to the absence of an analogue of 
the Boltzman-Gibbs approach and in spite of considerable recent progress [Q. 
Therefore non-equilibrium systems have to be specified by some defining dy- 
namical rules which are then analyzed. The topic has received a lot of attention 
and many reviews exist, e.g. 0-0 • 

One of the interesting and important examples of the non-equilibrium sys- 
tems is the one-dimensional reaction-diffusion processes, which have application 
in various fields of physics like study of the shocks noisy Burgers equation 
polymers in random media |10j. traffic models and biopolymerization |12j . 
As these systems are interacting systems with iV-particle, even simple models 
may pose a formidable problem if one wants to approach them analytically. See 
|13j-|16| for more recent references. 

The simplest reaction-diffusion process is the totally asymmetric simple ex- 
clusion process (TASEP). In this model, each lattice site is occupied by at 
most one particle and all particles can only hop with equal rate to their right- 
neighboring sites, if these sites are not occupied. TASEP has been studied in 
|17j by introducing a master equation which describes the evolution equation of 
the particles when they are not in neighboring sites, and a so-called boundary 
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condition, which specifies the situation in which the probabilities go outside the 
physical regions. This happens when some of the particles are in adjacent sites 
and the master equation can not be applied to them. It has been shown that the 
model is integrable in the sense that the iV-particle S'-matrix is factorized into a 
product of two-particle S-matrices. The coordinate Bethe ansatz has been used 
in this proof. 

The interesting observation is that if one chooses other boundary conditions, 
with the same master equation, one can in principle introduce other interactions 
(besides diffusion to right-neighboring sites), which may be integrable in the 
abovementioned sense. This is what is first done in ^B], m which the so-called 
drop-push model has been studied by this method. In this model the particle 
hops to the next right site, even it is occupied. It can hop by pushing all the 
neighboring particles to their next right sites, with a rate depending on the 
number of these particles. Some other generalization of TASEP can be found 



The generalization of one-species reaction-diffusion processes to p-species is 
an important task. The main problem in this generalization, besides introducing 
a set of suitable boundary conditions to model an interacting system, arises 
from the above mentioned factorization of ./V-particle scattering matrix. It was 
shown in |22| that in order that a more-than-one species system be solvable, in 
the sense of the Bethe ansatz, certain relations should be satisfied between the 
rates. These relations can be written as some kind of a spectral Yang-Baxter 
(SYB) equation. By this method, all the solvable two-species reaction-diffusion 
models, without annihilation and creation reactions and with equal reaction 
rates, have been obtained in |22j . 

The multi-species generalization of the reactions considered in [221 nas been 
studied in [2H], and the drop-push reaction of ^B] has been generalized to p- 
species in [23] . The most general totally asymmetric reaction-diffusion processes 
has been recently studied in |25|. These processes are 



where the dots indicate the other drop-push reactions with n-adjacent particles, 
in which in the meantime the types of the particles can also be changed. These 
latter reactions are called the extended drop-push processes. It has been shown 
that the reaction rates of processes Q must satisfy some specific constraints, 
in order that we have a set of consistent evolution equations. Also the corre- 
sponding two-particle S'-matrices must satisfy the SYB equation. Some classes 
of the solutions of these equations have been discussed in |25| . 

In all of the above studies, only the totally asymmetric exclusion processes 
have been considered, i.e. the particles can only diffuse to their next right neigh- 
boring sites. If one wants to consider the left and right diffusions simultaneously, 
one must consider a more general master equation with suitable boundary con- 
ditions and then seek the situations in which the model is integrable. In 17 , 
one-species model with only simple diffusion to left and right (i.e. partially 
asymmetric) has been considered, and in |26|. the one-species partially asym- 
metric drop-push model has been studied. Finally a two-species model in which 
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the particles, besides diffusion to left and right, have exchange-reaction has been 
studied in [2T| . 

In this paper we want to study the most general p-species integrable models 
with partially asymmetric reaction-diffusion processes, which all the previous 
studied models are the special cases of them. These general models may have 
some or all of the following reactions: 

A Q -> ®A a with rate D R , (2) 



%A a -> A a % with rate D Ll (3) 



A a A p -> A 7 A S with rate E"f?, (4) 



and 



A a Af3® -> %A^A S with rate R"%, 



$A a Af3 -> A 7 A 5 with rate L"f , 



(5) 



(6) 



In above equations a,fl, ■ ■ ■ = (1, • • • ,p), stands for vacancy, and dots in eqs.JSJ) 
and JJjJ indicate the drop-push of n-adjacent particles to right and left sites, 
respectively, in which in the meantime the types of the particles can also be 
changed. We call interactions JSJ and (JSJl as right-drop-pushing and left-drop- 
pushing, respectively. We show that there are three distinct models which are 
integrable and each of these models contains reactions (0 and and one or 
two of the reactions flj l to 

The sheme of the paper is as follows. There are two kinds of boundary 
conditions that can be generalized to p-species cases. In section 2, we generalize 
the first kind of boundary condition, which was introduced in |26) . to the most 
general p-species case. Using the law of conservation of number of particles, 
it is shown that there exists five constraints that must be satisfied by reaction 
rates of eqs.||2Jl-|©, m order to have a set of consistent evolution equations to 
express the interactions ©-©• But it is seen that there is no solution for these 
constraints. The situation does not change even if we relax one of the constraints 
by including the annihilation processes. Therefore one can not explain all the 
reactions 0-® by this method. But it will be shown that we can have two 
distinct models. In the first type model the reactions are eqs.(0, 0, an d 
(0) and in the second type the reactions are eqs.©, ©, 01 and ©. 

The second kind of boundary condition, which was used in ^7] and |2Z] , is 
generalized to the most general p-species case in section 3. We show that the 
resulting consistent boundary condition can explain the reactions (01, 10 and 
(@J. This is the type 3 model. It must be mentioned that the type 3 model 
is not a subclass of types 1 and 2 and is a new distinct one. In section 4 we 
investigate the Bethe ansatz solution for these models and discuss the solutions 
of the corresponding SYB equations. We see that the 5-matrix of type 3 model 
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is much more involved than two other ones and therefore only some special 
classes of solutions of its SYB equation can be obtained. Finally we study the 
conditional probabilities of these models and in special two-particle sector, we 
obtain the exact expressions. 

2 First kind generalization 

Consider ap-species system with particles A%, A2, ■ ■ ■ , A p . The basic objects we 
are interested in are the probabilities P ai -a N ■ • • , %n] t) for finding at time 
t the particle of type a.\ at site x\, particle of type 012 at site X2, etc.. We take 
the physical region of coordinates as x\ < X2 < ■■■ < xn- The master equation 
for a partially asymmetric exclusion process is 

N 

— 1, x i+ i, ■ ■ ■ , xn; t) 

1=1 

JV 

+ 1, Xj+i, • • • , xn', t) 

i=l 

-NP otl ,.. olN (x 1 , ■ ■ ■ ,x N ;t). (7) 

This equation describes a collection of N particles, diffusing to the next-right 
sites by rate Dr and to the next-left sites by rate D^. In eq.JZJ) we have used 
a time scale so that 

D R + D L = 1. (8) 

This master equation is only valid for Xi < a;i+i — 1. For xi = Xi + \ — 1, there will 
be some terms with Xi = 2^+1 in the right hand side of eq.Q which are out of the 
physical region. But one can assume that Q is valid for all the physical region 
xi < Xi+i by imposing certain boundary conditions for Xi — x^+i. Different 
boundary condition introduces different interactions for particles. Following the 
argument which have been given in it can be easily seen that the master 
equation JJJ) leads to following relation for two-particle probabilities: 

y^}D R P aia2 (x,x; t) + D L P aia2 {x + l,x+ l;t)] ~^P ai0t2 (x,x + l;t) 

X X 

^ ] Pu\Qf2 {x,x,t) — S ' P ai a 2 {x, X + 1, t). (9) 

X X 

This equation leads us to take P ai a 2 ( x , x; t) as linear combination of P^fo ( x , x+ 
l;t) and P / 3 1( g 2 (x — l,a;;t)'s as the only choice for having a consistent set of 
evolution equations in more-than- two-particle sectors |25j . Therefore the most 
general boundary condition is 

P aia2 (x,x) =^^ 2 2 P ft/32 (x-l,a ; ) + ^ C ^ 2 P /3l/32 (a ; ,x + l). (10) 

(3 

(3 stands for (fiifo) and b and c are p 2 xp 2 matrices determine the interactions. In 

the probabilities appear in eq. l(TU|) , we have suppressed all the other coordinates 

and the time t for simplicity. In fact P ai a 2 ( x i x ) ■— P-yi---y i a 1 a 2 'y i+3 ---y N {xi, ■ ■ ■ , Xi, x, x, Xi+3, ■ ■ ■ ,xn)- 

In the first step, let us exclude the creation and annihilation processes (it can 



dt a 



,(x\, ■ ■ ■ ,x N ]t) 



dt 



Paiai^l^^ 1 ) = 



x 2 x±<x 2 
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be shown that one can not study the creation processes by this method, so in 
fact in this step, we exclude the annihilation processes). Since the number of 
particles is constant in time, summing over ai and a-i makes the left-hand side 
of © zero and results: 

-EE p ^ (*> x + !) + E E ( E( 6 + ) x+i) = o, (ii) 

in which ea. (|10|l has been used. Clearly ea. (|llfl gives: 

E( 6 + C )«i«2 = 1 constraint (I). (12) 

a 

Note that in p = 1, the boundary condition (|10|l and constraint l|12|) reduce to 
those considered in [2SI - Also in the case of totally asymmetric processes in 
which Dl = 0, our problem reduces to one considered in |25j . Following the 
same steps as we first consider P aiCt2 (x,x + 1). Using eqs.0 and (ffTTTf) - it 
is found 

P Q1Q3 (x,x + 1) = D^ aia2 (x-l,x+l) + -D£P aiQ3 (x,x + 2) 

+ D R E W^ft " !> *) + ^ E ^.^Aft + 1, i + 2) 

^ C ^ 2 2 + ^ 6^ 2 2 )P PlPa (x, x + 1) - 2P Q1Q2 (x, x + 1) 

= D R P aia2 (x - l,x+ 1) + D L P aia2 (x, x + 2) + D R ^bi\^ 2 P M2 {x - 1, x) 

P 

+Dl 53 Pfcfl, (x + 1, x + 2) + 53 (Dh^* + JPftfc (*> s + l) 

-p* + £>l + E(^ c ft& + ^TO + E(^X + ^« 2 )]p Qia2 (x,x + 1), 

(13) 

in which we use eqs.ljHJ) and (|12|l . The latter can be written as: 

ftlft 2 _ 1 _ \ ^ \ ^ „«1«2 /I y|\ 

/3 /3#a 



or 



Lft 1 ft 2 _ 1 \ ~" 1.0102 aift2 /IC-l 

°ftift 2 - 1 Zv °/3i/3 2 ~ Zv C ft/32 • ^ i0 J 

It is seen that the evolution equation ifHfll describes the following two-particle 
interactions: 



(16) 





-> 0A Q 




-+ A Q 


A a Ap - 


-> A^Ag 


A a A fj (D - 


-> %A 1 A S 


®A a A (j - 


■* A y A s ® 
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To study the consistency of our formalism and also deriving the more-than-two 
particle interactions, we considerP Q , 1 ... Qri (x, x + 1, • • • , x + n — 1). In n = 3, we 
encounter two boundary terms P aia2 a 3 (x,x + l,x + 1) and P ai a 2 a 3 {x + 1, x + 
1,35 + 2). Using the first one becomes: 

Pa 1 u 2 a 3 {x, X + 1, X + 1) = ^7172/83 (^ — 1, X + 1) + ^1^2-^7172/33 (Xj ZE + lj 21 + 1)] 

+ E C £«3 P « 1 ^3^^ + 1 ^ + 2 ) ( 17 ) 

which describes the boundary term P QlQ2Q3 (a;, x + 1, x + 1) as a linear combina- 
tion of other boundary terms, i.e. Pf^fh (x, x + 1, x + l)'s. As has been shown 
in |25], the only consistent solution to this problem is the vanishing of these 
terms in the right-hand side of ea . Ij 1 7|1 . which results: 

E^fttt=° constraint (II), (18) 

ft 

or 

(1 ®&)(c<8) 1) = 0, (19) 
in which 1 stands for the p x p identity matrix. The second boundary term is 

P aia2a3 ( x + l,x + l,x + 2) = J2bP\i 2 2 Pp ll 3 2a3 (x,x + l,x + 2) 
_^ 
+ E c tt HIZ P ^ 3 (x + l,x + l,x + 2) + c}ll\P 0m (x + 1, * + 2, z + 3)], (20) 

01 

which again leads us to take 

E c S* & ft£=° constraint (III), (21) 

ft 

or 

(c(g>l)(l<g>&) = 0. (22) 

Assuming constraints (fT%|) and (f^Hl and using eqs.Q and (fTO)) . P,s(a;, jc+l, x+2) 
is 

P s {x,x + l,x + 2) = D fl F 5 (x-l,x + l,3j + 2) + L> L P 5 (x,a; + l,x + 3) 

E tt^ft/W* - 1, x + 2) + £>i ft^ftA (*> * + 2, z + 3) 

+ E (^i* + (i, a; + 1, a: + 2) 

+ ^ (I> fl c** + D L bif a %)P aM (x,x + l,x + 2) 

+ D R E 4^ - 1. a: + 1) + D L E 4 P 7^ + 1, z + 2, z + 3) 

-\Dn +D L + D R £ 6££ + Z? L E^ft + ^ E & ftft + ^ E c ftft 

/5 |9 /3 

+ E ( D « c ftft + ^ftft) + E (^ c ftft + D L b;^)]P s (x, x + 1, x + 2), 



6 



(23) 

in which we have used eas. p4(l and ljl5|l for diagonal elements of matrix D R c + 
Dib. bZ and cl arc defined as following 

4 = E b l 1 X 3 ^ ( 24 ) 

7 

r l _ r lxP „7273 fort 
7 

Looking at source terms of ea. (|23|l . it is obvious that they describe the reactions 
Ultjfl and the following three-particle drop-push reactions: 

A 71 A 72 A 73 —> $A ai A a2 A a3 with rate D R bZ, (26) 

and 

0A 7l A- y2 A 73 — > A ai A a . 2 A a . 3 % with rate D L ct. (27) 
The sink terms are consistent with this description, provided 

b txpl constraint ( IV ), ( 28 ) 

P 

^Z° a plZ constraint (V). (29) 

P 

By calculating other Ps(x, x + 1, • • • , x + n — l)'s it can be shown that we need 
not any more constraints and therefore the master equation Q with boundary 
condition (|Tfj|) and five constraints (I)-(V) can consistently describe the following 



reactions: 

A a % — ► 9A a with rate D R (30) 

$A a — y A a 9 with rate D L (31) 

A a A ^A 7 A s with rate D R c$ + D L b$ (32) 

A ao ---A a J^%A^ ---A ln with rate Dr(&„_i,„ • • •6o,i)^.'.'.'^ , ) (33) 

0A« O • • ■ A an -> A 7o ■ ■ • A 7 „0 with rate r»i(co,i • ■ • Cn-i.n)^"^". (34) 

In above equations we use the following definition for bi^k+i and Ck,k+i '■ 

a k ,k+i = 1 ® • • • ® 1 ® ®1 <g> • • • ® 1. (35) 
fc,fc+i 



O 0~2^°017 O P203 



01 



and 



E„5 _ «i7 „«2«; 

C /3 ~ 2-^, L 0102 103 
01 
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Note that for Dl = 0, the five classes of the above reactions reduce to three 
ones discussed in In [35], the constraints between reaction rates are three 
relations (I), (II), and (IV). Note that at Dl = 0, the constraints (III) and (V) 
do not appear since the multiplication factors of their corresponding terms in 
evolution equation is Dl , which is zero. 

To find the set of solutions of five constraints (I) to (V), one can consider 
the solutions of equations (I), (II) and (IV), that is the solutions derived in [25] . 
and then considers the subset of them satisfys (III) and (V) . We must note that 
in our models, the diagonal elements of matrix c are the reaction rates of the 
last line of ea. (|16f) and must be positive. This is in contrast to the case studied 
in |25) in which the diagonal elements of c can be negative. 

We can also follow another approach. That is trying to find the solution of 
equations (II) to (IV) and then seek ones which satisfy relation (I). As these 
relations are rather complex, we can not completely solve them for arbitrary p, 
but we try them as much as possible. 

As all the matrix elements of matrices b and c are reaction rates, they can 
not be negative, so the only solution of ea. l|18|) is: 

<&tt=° (without sum over /3 2 ). (36) 

This relation has two following solutions ( for each /3 2 ) 

cl\%=0 and 6^ 3 =0. (37) 

So for each f} 2 we have two solutions, and as /3 2 runs from 1 to p, we have 2 P — 2 
set of solutions for constraint (II). We exclude two of the solutions in which 
all of the elements of c or b is zero , since we look for the situations in which 
6^0 and c ^ 0. We will later study the cases b = or c = in which the 
number of independent classes of reactions (j30|l - l|34(l reduces to four. By the 
same argument, the solutions of ea. l|2T|l are (for each /? 2 ) 

c**=0 and &££=(), (38) 

and therefore we again have 2 P — 2 set of solutions for constraint (III). Note 
that from 2 P — 2 solutions of constraints II ( and III ) only (p — 1) of them are 
independent, that is does not transform to each other under interchanging of the 
labels of the species of the particles. So the number of independent solutions 
of constraints (II) and (III) are (p — 1)(2 P — 2). For example in p = 2, the 
independent solutions of (II) and (III) are 

{ ci\r = o, bi% 3 = o, ci\\ 2 = o, bf:: = o>, 

{ cl\f = 0, bfx = 0, 4\l 2 = 0, bj™ = 0}, (39) 
which can be written as 
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respectively. We label the states as |1 >= (1,1) , 2 >= (1,2) , |3 >= 
(2,1) and |4 >= (2,2). Putting eqs.gUJ and (gTJ into the constraint (IV) 
and (V) (eqs.lJH and J2SJ) results 



6 = 



and 
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(42) 



(43) 



respectively. Now eq. (|12f) (constraint I) says that the sum of the elements of each 
column of matrix (b + c) must be one, which unfortunately does not satisfy by 
(|42|l and 143JI . The sum of the elements of the second column of (b + c) of ea. (|42|) 
and the first column of eq. l|43ll are greater than (or equal to) 2. So reactions 
(I30|l - H34ll have not any representation in p = 2, the situation which we expect 
to be true for other p's. For example in p = 3, constraints (II) and (III) have 12 
independent solutions, which are in two categories: the number of constraints 
on 6/ and c fc ''s are equal (three on 6/ and three on c fe '), and one which these 
numbers differ, i.e. 4 and 2. As an example of the first category, we consider 



the case in which c„ 



0,c 



a 2 



0,c 



/31 



0,6 



3 7 



Oh k 
u i u la 







h k 
U 3a 



0. It 



means that in matrix c, the rows 1,2,4,5,7,8 and columns 2,5, and 8 are zero, so 
it has 18 non-zero elements, and b is a matrix in which the rows 1,2,3,7,8,9 and 
columns 7,8, and 9 are zero so it also has 18 non-zero elements. Putting these 
b and c matrices in constraint (IV) and (V) results two following solutions for 
each b and c: 



Cl 


{C31 


= C34 


= C37 


= l,c 33 


and 


C36 


arbitrary}, 


C2 


{C93 


= c 96 


= C99 


= l,c 9 i 


and 


C94 


arbitrary}, 


h 


{^41 


= 642 


= 643 


= 1,644 


and 


645 


arbitrary}, 


b 2 


{&54 


= &55 


= &56 


= 1,651 


and 


652 


arbitrary}, 



(44) 

in which we only write down the non-zero elements. It can be easily seen that 
none of the combinations b\ + c\, b\ + C2, 62 + Ci, 62 + C2 are acceptable in 
the sense of constraint (I), as at least the sum of the elements of one of the 
columns of these matrices are greater than (or equal to) 2. We have checked 
that the same situation arises in other 11 solutions. So again in p — 3, we have 
no representation. We can not generally prove this, but we believe that the set 
of constraints (I)-(V) have no solution for arbitrary p. 

One may suppose that if we somehow change the constraint (I), then it may 
be possible to find some solution for our equations. So we add the annihilation 
processes to our previous interactions. Note that these interactions appear only 
in the sink terms of the evolution equation, as if we consider the initial state with 
n particles, no annihilation processes can lead to a n-particle state at any other 
time t. So if we change the constraint (I) to: (as we have not the conservation 
of particles) 

E( 6 + C tt =1 - W (45) 
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and using it in calculation of P aiCC2 (x,x + 1), we find the same equation as 
(|13|l . except an extra term X aia2 P aia2 (x 1 x + 1) which is added to sink terms. 
So X ai a 2 is the sum of the rates of all annihilation processes with initial state 
(aia^) and therefore is a positive quantity. Therefore adding the annihilation 
processes to interactions (|rt()l) - H34|) means that the sum of the elements of each 
column of (6 + c) can now be less than or equal to one. But as we have shown in 
eas. <|42[) . (|43|) and l|44|) . the sum of the elements of some of the columns of (b + c) 
in these examples are at least 2, which differs from what is suggested by ea. (|45ll . 
In brief, including the annihilation processes can not alter our result and the 
set of processes (I3C)I> - (|34l> have no representation, with or without adding the 
annihilation processes. 

Now it is interesting to note that even if one of the matrices b or c be equal 
to zero, we have yet all four desired reactions: Diffusion to left and right, two- 
particle reactions A a Ap — > A^Ag, and the extended drop-push reactions, which 
the latter occur only in one side (left or right). These are almost the general 
reactions that one can study in this framework. Let us check the constraints in 
these cases. 



2.1 Type 1 model 

Take c = 0. Eq.(12) becomes 

Ett = L ( 46 ) 

a 

Constraints (II), (III) and (V) are satisfied trivially and constraint (IV) is also 
satisfied: using eq. l|46|l . both sides of (IV) become one. Therefore master equa- 
tion J7J with boundary condition 

P aia2 {x,x) =]>>f 1 1 * 2 P /3l/32 (x-l,x), (47) 



and constraint (|46J) . describe consistently the following reactions: 

A a % -> %A a with rate D Rl 

$A a -» A a % with rate D L , 

A a Ap — > A y As with rate Dlo^, 

A ao ---A a J -» 0A 7o ---A 7 „ with rate D R (b n - 1<n ■ ■ ■ b ,i)%:::%-{^) 

2.2 Type 2 model 

In the same way, for b = it can be seen that the master equation (0 with 
boundary condition 

P aia2 (x,x) ^cf*%( IjI + 1), (49) 


and constraint 

E c ^ = 1 > ( 5 °) 

a 

describe successfully the reactions: 

A a % -> 0A Q with rate D R , 
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%A a -> A a 9 with rate D L , 
A a Ap -> A 7 A 5 with rate D R c*g, 
%A ao ---A an -» A 7o ---A 7 „0 with rate D £ (co,i • • • c n _i, n )^.7.^. (51) 

The condition of solvability of these models will be discussed in next sections. 

3 Second kind generalization 

By noting the first line of eq.JHJ), it is seen that ea. (|l()Jl is not the only possible 
p-species boundary condition. In fact, one can instead consider the following 
boundary condition : 

D R P aia2 (x, x)+D L P aia2 (x+1, x+l) = ift&Pfrfh (x-1, C ' 3 S\ P P^ ( x > x + 1 )- 



(52) 

This is the multi-species generalization of the boundary condition considered in 
QU and [2Zj. 

To study the interactions introduced by Q and \\W1\\ , we must again consider 
Pat--a n (x,x + 1, • • • , x + n — 1). In n = 3, we encounter the boundary term 
D R P ai a 2 a 3 (x, x + 1, x + 1) + DLP aia20l3 (x, x + 2, x + 2), where using l|52|) results 
E/3 b1f a 3 3 P ai f3 2 i3 a (x, x, x+l)+Y,p c ala 3 p <*i0203 i x , x +l, x+2). But the first term 
P ai /3 2 p 3 (x, x, x + 1) can not be written in terms of physical probabilities, since 
in this case only the linear combination D R P aiCt2 ...(x,x, • ■ • ) + D] J P aia2 ...(x + 
l,x + 1, • • •) can be written in terms of physical function fea. l|52|l 'l. This is in 
contrast with the case studied in section 2. The only solution to this problem 
is taking 

6 = 0. (53) 

So our second kind p-species model is defined through the master equation l|7| 
and the following boundary condition : 

D R P aia2 (x, x) + D L P aia2 {x + 1, x + 1) = 4\i\P0x0 2 (x, x + l). (54) 



Conservation of number of particle gives 

a 

and calculating P aia2 (x, x + 1) results 

P aia2 {x, x + l) = D R P aia2 {x - l,x+ 1) + D L P aia2 (x, x + 2) + (&[a 2 P0i02(x, x + l) 

0+a 

~{D R + D l + Y, c a ^ 2 2 )P QlQ2 (x, x + 1). (56) 

0^a 

This equation describes the following reactions as source and sink terms 

A a 9 -» $A a with rate D R , 
%A a -» A a % with rate D L , 
A a Af3 -> A 7 A S with rate c°f & . (57) 

Calculating other P ai ... arl (x, x + 1, • • • , x + n — l)'s confirms these reactions 
without any further constraint. So the type 3 model is defined by master 
equation J7J|, boundary condition (|54|l . constraint H55[) . and reactions (|57|l . 
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4 Bethe ansatz solution 

Until now, we have constructed a consistent formalism to study some reaction- 
diffusion processes. Now we want to solve the resulting evolution equations and 
check the solvability of these models. To solve the master equation Q, we 
consider the following Bethe ansatz 

P ai -* N (x; t) = e~ ENt ^ ai ... aN (x), (58) 

with 

#(x) =^A ff e 4CT ( p K (59) 

a 

^ is a tensor of rank N with components ip ai ... aN (x) and the summation runs 
over the elements of the permutation group of N objects 28, 29 . Inserting i|58|) 
in J2J), results 

N 

E N = J2( 1 - D R e ~ iPk -D L e ipk ). (60) 
fe=i 

Inserting l|58(l in boundary condition 147J) gives 

*(• • • ,x k = x,x k +i = x, • • •) = /, . • • ,x k = x - l,x k +i = x,- ■■), (61) 
which using ljS*9"|) results 

[1 _ e — M bk fc+1 ] Aff + [1 - e - ia ^+^b ktk+1 ]A aak = 0. (62) 
a k is an element of permutation group which only interchanges p k and Pfc+i: 

a k : (pi, • ■ • ,p k ,p k+1 , ■■■ ,p N ) -> (pi, • ■ • ,p k+1 ,p k , ■ ■ ■ ,p N ). (63) 
Ea. (|62|l gives A ja|! in terms of A CT as following: 

= 4!fc+i( CT (^)^(Pfe+i))A CT , (64) 

where 

S^k+ifci,^) = 1® ■••<& 1®S w (z 1 ,z 2 )®1®---®1, (65) 

fc,fc+i 

and S'^^zi, z^) is the following p 2 x p 2 matrix 

S (1) (^ 2 ) = -{l-z^b)-\l-z^b), (66) 

in which 2^ = e lpk . The same procedure for boundary conditions (|49|l and (|54J) . 
i.e. the type 2 and type 3 models, results 

S^( Zl ,z 2 ) = -(1 - zic)~\l - z 2 c), (67) 

and 

S (3) ( Zl ,z 2 ) = -(D R + Zl z 2 D L - z 1 c)- 1 (D R + Zl z 2 D L - z 2 c), (68) 

respectively. Ea. (|64|l allows one to compute all A CT 's in terms of Ai (which is 
set to unity). 
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As the generators of permutation group satisfy Ok(Jk+iOk 
one also needs 

A = A 

JrL <^k<^k + l<^k ^fk + l&k&k + l' 

This, in terms of S'-matrices becomes 

£12(22, 23)523(21, z 3 )Si 2 (z 1 , z 2 ) = £23(21, 22)512(21, 23)523(22, 23). (70) 
In the terms of i?-matrix defined through 

5fc,fc+l =: Iik,k+lRk,k+li (71) 
where II is the permutation matrix, ea. (|7()|) is transformed to 

#23(22, z 3 )R 13 (zi, z 3 )R 12 (zi,z 2 ) = Ru(zi, z 2 )Ri 3 (zi, z 3 )R 23 (z 2l z 3 ). (72) 

This is the spectral Yang-Baxter equation. 

The Bethe ansatz solution exists, if the scattering matrix satisfies Q7()|l. In 
other words, the matrix b in l|66[l and c in (|67|l and (|68[1 is acceptable only if the 
resulting S*-matrices satisfy (|70|) . This is a very restricted condition and needed 
for having the solvability. 

The S'-matrices (166(1 and (16711 are exactly the ones considered in j23] and , 
respectively. Using the fact that is a binomial of degree one with respect to 
z± 1 = e~ lpi and is of degree one with respect to z 2 , it can be shown that 
SYB equation (|7UJ) for and reduces to 

b 23 [b 23 ,b 12 ] = [b 2 3,bi 2 ]bi2, (73) 

and 

C12[C12, C23] = [C12, C 2 3]c 2 3, (74) 

respectively [2B1 E] . Note that although the above equations are much simpler 
than ea. (|70|l . but they are very complicated yet. In p-species, each one is an 
equality between two p 3 x p 3 matrices which results a system of p 6 equations to 
be solved for p 4 — p 2 elements of b (or c) , which may or may not have solution 
fca. (|46|l and H50[l reduce the number of independent elements of b and c to 
p 4 — p 2 ). The general properties of the solutions of eq.JJSJ and ea. (!74|l have 
been discussed in [21] and [2Sj , respectively, which can be directly used here. In 
other words, for every solution of ea. l|74|l . there exists a corresponding solvable 
model which have been discussed in [231 , i- e - 

A a % -» %A a with rate 1, 
A a Af3 -» A 7 A S with rate c"f . (75) 

and a type 2 model with reactions written in ea. l|51|) (note that the reactions 
(|75[1 are a subset of 1)51(1 with Dl — 0). The same is true for solutions of l(73|l . 
They can describe the following solvable model (discussed in |24|): 

A a 9 — > ®A a with rate 1, 
A ao ---A a J -► %A 1Q ---A ln with rate (6„_i,„ • • • 6 0> i)" ° (76) 

and a type 1 model with reaction l|48(l (again at Dl — 0, (|48|1 reduces to i|76|l ). 

The reasoning which leads the SYB equations of S^ 1 ' and to l|73|) and 
(I74|l does not work for since it is not a binomial of degree one with respect 



= (Jfe+lCfcCfc+l, so 

(69) 
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to z\ or Z2, in fact it contains all powers of Z\ and z^. So obtaining the solutions 
of J701 for S ( V is more difficult than for and S^ 2 ), even in the simplest 
case p = 2. In p = 2 case we encounter a system of 64 equations that must be 
solved for 12 non-diagonal elements of c (the diagonal elements are determined 
by ea. 155(0 . The solution must be momentum-independent (independent of z\, 
zi and Z3) and non- negative. We can not solve this equations generally ( taking 
all Cij ^ ) by standard mathematical softwares and therefore restrict ourselves 
to some specific cases. For example taking 



or 



Cll 












C14 
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C21 




1 







C24 




en - 


C21 
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- C14 - 


C24 


















/ 
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C21 




1 







C24 




C31 







1 




C34 
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C31 
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- C24 - 


C34 / 



(77) 



(78) 



V 1 

which are the four-parameters cases, one obtains two solutions 



( ° 

Dr 



\ 

1 D L 
1 D R 



(79) 



or 
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^Dr 
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(80) 



(81) 



as some three-parameters cases, we find four solutions 



/I 





V 



1 - C42 1 - C43 






C42 






C43 






1 / 



where each of C42 and C43 are either Dr or Dl. Taking A\ 
the interactions introduced by lf79)l. for instance, are 



(82) 



A and Ao = B, 



A® 


Dr 


0-4 


Bfb 


Dr 


%B 


9A 


Dl 


Ail) 


0B 


Dl 


BID 
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AA ^ AB 

BB °i AB 

AA °i BA 

BB 5$ BA. (83) 

The model built on the reactions l|83|l is integrable. 

Assuming that the solvability condition (|7[)|) is satisfied, it is easy to see that 
the conditional probability (the propagator) is 

f rl N r> 

C/(x; t\y; 0) = J ^-^e^e^ £ A^^, (84) 

where the integration region for each pi is [0, 2tt] and A\ = 1. The singularity 
in A a is removed by setting pj — > pj + ie, where one should consider the limit 
e — > + . Using this propagator, one can write the probability at the time t in 
terms of the initial value of probability: 

|P(x;*))=^*7(x;i|y;0)|P(y;0)}. (85) 
v 

Note that although and are similar to ones considered in [21] and |23| . 
but the propagators f/Wand are different since the energy spectrum of our 
models differs from those considered there. In Dl — 0, our results must coincide 
with those obtained in 123 E3] ■ 

For the two-particle sector, there is only one matrix in the expression of 
l/«s {b in [/« and c in and U^). So it can be treat as a c-number. Using 
calculation similar to what has been done in |23|-|25|. one arrives at: 

°° 7-m T~ l xi-yi+n,x 1 -y 1 +2n nmn :I i"9!+" 1 ii!-V2+2m 

' n - yi + n)! m!(x 2 - y2 + raj' 

n,m— v 7 v 7 

°° f) X2 ~ 2/i+ n +a?2— 2/i+2ti jjm jj x i—y2—l+rn, X i—y2—l-\-2rn 
j \. L R L i?, 



^ n!(x 2 - ?/i + n)! m!(:ci — 1/2 — ' + "i)! 

+ (86) 



and 



^)(, 1 ,, 2 ;% 1 , 2/2 ;0) = e- 2 * £ V v 

' n!to — yi + n)! m!(x 2 - y 2 + raj! 

n,m— v y 

£jj n!(a; 2 - ?/i + / + ra)! m!(xi - yi + m)\ 

xc l [-l+ ^ -c]}. (87) 

aci -y 2 + m + l u y ' 

Similarly one can obtain a more lengthy expression for . Note that at 
D L = 0, eqs.JEHJ and JH3 lead eqs.® of [22 and (JHDJl of 23 . respectively. 

To investigate the large-time behaviours of the probabilities £/W, U^- 2 \ and 
L/C 3 ) it is useful to decompose the vector spaces on which b (in type 1 model) 
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and c (in types 2 and 3 models) act, in two subspaces invariant under the action 
of b(c): the first subspace corresponding to eigenvalues with modulus one, and 
another invariant subspace. For types 1 and 2 models with conditions (|46() and 
1)50(1. as all the elements of matrix b(c) are non- negative, the second subspace 
corresponds to eigenvalues with modulus less than one. By focusing on type 1 
model, this decomposition can be done by introducing two projectors Q and i?, 
satisfying 

Q + R = 1, 
QR = RQ = 0, 
[b, Q] = [b, R] = 0. (88) 

Q projects on the first subspace and R projects on the second. Following |23| . 
we multiply C/ (1) by Q + R := 1: 

U {1) (x; i|y; 0) = U {1) Q + U {1) R. (89) 

In the terms multiplied by R, one can treat b as a number with modulus different 
from one. So the integrand in l|84|l is non-singular at points pj = 0, which have 
the main contributions at large times. Putting pj = 0, we have « -1 and 
A a w (-1) M , and eq.(JHH) results 

the second term of E/« = — { e -{[*i-vi-(D R -D L )tf+[x 2 -y 2 -(D R -D L )tf}/(2t) 

_ e -{[^l-y2-(DH--D I ,)t] 2 + [ 2;2 -y 1 -(D„-I5 I/ )t] 2 }/(2t)| jR 



which is independent of 6. So at large time, the second term of tends to 
zero faster than t~ l and the leading term in U^ 1 ' , which is order t -1 , does not 
involve the second term. 

If the only eigenvalue of b with modulus 1 is 1, then bQ = Q and has a 
simple behaviour at t — » oo: 

uW( Xl , X2] t\ yi ,y 2 ;0) = e~ 2t ^ — — — 

' nUXi - yi + n)! m!(a;2 - ife + my. 

n,m—0 v ' 

00 ^ j~)n jjX2-yi+ri£ X2 -yi+2n j-jm j-)Xl-y2-l+m_f_x 1 -y 2 -l+2ni 

^ n!(x 2 - ?/i + n)! m!(a;i - j/2 - I + m)\ 

x[-l + ^^]}Q. (91) 

This is simply the propagator of a single-species model with diffusions to right 
and left and drop-push to right (i.e. the A = case of the reactions studied in 
[23 and [H]), multiplied by Q. In fact eq.JS2J is A = case of eq.(30) of gEj, 
times Q. 

For U {2 \ the same decomposition leads the ea. (l9UII for its second term and 
in the case cQ — Q, tends to (|87l) . with c = 1, times Q, at t — ■> oo. The 
resulting one-species model is /j. = case of the reactions studied in and 
[S3j. For t/( 3 \ we again find <HI3J) and the one-species partially asymmetric 
simple exclusion process of |17| . 
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